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We discuss the phenomenology of the lightest stop in models where R-parity is broken by 
bilinear superpotential terms. In this class of models we consider scenarios where the R-parity 



> 
O 

breaking two-body decay t\ —* t + b competes with the leading three-body decays such as t\ — > 

""xt" \ , n 

W + bxi- We demonstrate that the R-parity violating decay can be sizable and in some parts of 
the parameter space even the dominant one. Moreover we discuss the expectations for t% — > fi + b 
and t\ — > e + b. The recent results from solar and atmospheric neutrinos suggest that these are 
as important as the T + b mode. The t\ — > / + b decays are of particular interest for hadron 
colliders, as they may allow a full mass reconstruction of the lighter stop. Moreover these decay 
^ ■ modes allow cross checks on the neutrino mixing angle involved in the solar neutrino puzzle 
complementary to those possible using neutralino decays. For the so-called small mixing angle 
or SMA solution t\ — > e + b should be negligible, while for the large mixing angle type solutions 
all t\ — ► / + b decays should have comparable magnitude. 
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I. INTRODUCTION 



The search for supersymmetry (SUSY) [|T], |2[ plays an important role in the experi- 
mental program of present and at future colliders, e.g. the Tevatron, LHC, or an e + e~ 
linear collider. Therefore many phenomenological studies have been carried out in re- 
cent years (see e.g. |3|, f|, [5], |6| and references therein) focusing mainly on the minimal 
supersymmetric standard model (MSSM) 0. However, neither gauge invariance nor 
supersymmetry require the conservation of R-parity. Indeed, there is considerable the- 
oretical and phenomenological interest in studying possible implications of alternative 
scenarios || in which R-parity is broken || [K| [11], |12| . These theories are of particular 
interest as they lead to a pattern of neutrino masses and mixing angles |13| which can 
account for the observed anomalies in solar and atmospheric neutrinos [Q. In general 
the violation of R-parity could arise explicitly |15[ as a residual effect of some larger uni- 
fied theory [|l(J , or spontaneously, through nonzero vacuum expectation values (vev's) for 
scalar neutrinos || [11], [12], [Uj . In realistic spontaneous R-parity breaking models there 
is an SU(2) (g) U(l) singlet sneutrino vacuum expectation value (vev) characterizing the 
scale of R-parity violation [IT], [TE|, ID], ^U] which is expected to be in the order of 1 TeV. 



There are two generic cases of spontaneous R-parity breaking models to consider. In the 
absence of any additional gauge symmetry, these models lead to the existence of a physical 
massless Nambu-Goldstone boson, called majoron (J) which is the lightest SUSY particle, 
massless and therefore stable. As in the standard case in R-parity breaking models the 
lightest SUSY particle (LSP) is in general a neutralino. However, it now decays mostly into 
visible states, therefore diluting the missing momentum signal and bringing in increased 
multiplicity events which arise mainly from three-body decays such as x? —* ffv, where 
/ denotes a charged fermion [|!J. If lepton number is part of the gauge symmetry and R- 
parity is spontaneously broken then there is an additional gauge boson which gets mass 
via the Higgs mechanism, and there is no physical Goldstone boson In this case 

R-parity violating effects relevant for collider physics are conveniently parameterized by 
adding bilinear terms ejLjif 2 to the MSSM superpotential and corresponding terms for 
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the soft SUSY breaking part of the Lagrangian. Bilinear R-parity violation may also be 
assumed ab initio as the fundamental theory. For example, it may be the only violation 
permitted by higher Abelian flavour symmetries pi|]. 

Owing to the large top Yukawa coupling the stops have a quite different phenomenology 
compared to those of the first two generations of up-type squarks (see e.g. |22], ^3|] and 
references therein). The large Yukawa coupling implies a large mixing between ti and tn 
Pfl and large couplings to the higgsino components of neutralinos and charginos. The 
large top quark mass also implies the existence of scenarios where all MSSM two-body 
decay modes of t\ are kinematically forbidden at the tree-level (e.g. ti — > t Xi, bxt, tg)- 



In such case higher order decays of ti become relevant p5L EBL |27|, pOl: ii — > cx° 
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ti — > W + bx\, ti — > H + bxi, ^i — > blfb>i, ti — > buil + , where I denotes e,/i, r. Also 
4-body decays may become important if the 3-body decays are kinematically forbidden 
[ p8|j . In ^7\, |3(| it has been shown that in the MSSM the three-body decay modes 



are in general much more important than the two body decay mode. Recently it has been 
demonstrated that not only LSP decays but also the light stop can be a good candidate 
for observing R-parity violation, even if its magnitude is as small as indicated by the 
solutions to the present neutrino anomalies f| fTj], |32|, |33| . It has been demonstrated that 



there exists a large parameter region where the R-parity violating decay 

ti -> br 

is much more important than the conserving decays 

*i C Xl,2 

in scenarios where two-body decay modes are possible. It is therefore natural to ask if 
there exist scenarios where the decay t\ — ► b r is as important as the three-body decays. 
Note that in the R-parity violating models under consideration the neutral (charged) 
Higgs-bosons mix with the neutral (charged) sleptons [33]]. These states are denoted 
by S®, Pj, and for the neutral scalars, pseudoscalars and charged scalars, respectively. 
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Therefore in the R-parity violating case one has the following three-body decay modes: 

ti -> S£bu 3 

fx - 6 5? r+, 

ti -> 6P?r + 

ii — > 6 z>i Z + (7 = e,fj) . 

We will show that there exist regions in parameter space where ii — > b r + is sizeable and 
even the most important decay mode. In particular we will consider a mass range for the 
light stop t\, where it is difficult for the LHC to discover it in the MSSM due to the large 
top background . In contrast to the existing LSP decay studies in $, p models |]37|, |38 



which are mainly sensitive to the atmospheric neutrino anomaly parameters, the stop 
decay processes considered here are very sensitive to the solar neutrino parameters and 
therefore gives valuable complementary information. 

The paper is organized in the following way: in the next section we will introduce 
the model. In Sect. |II| numerical results for stop decays are presented. We first explore 
the extent to which the decay t\ — > br can be sizeable when compared with the 3-body 
decay modes. Moreover, we discuss the connections between the decay modes t± — > bl + 
and neutrino physics, in particular we discuss a possible test of the solution to the solar 
neutrino puzzle. In Sect. [IV] we present our conclusions. The appendixes contain complete 
formulas for the total widths of the three-body decay modes as well as for the couplings. 

II. THE MODEL 

The supersymmetric Lagrangian is specified by the superpotential W given by 



W = e ab 



2 ■ 
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where i,j = 1,2,3 are generation indices, a,b — 1,2 are SU(2) indices, and e is a com- 
pletely antisymmetric 2x2 matrix, with £12 = 1. The symbol "hat" over each letter 
indicates a superfield, with Q i: Lj, H 1: and H 2 being SU(2) doublets with hypercharges 
1/3, —1, —1, and 1 respectively, and U, D, and R being £77(2) singlets with hypercharges 
— |, |, and 2 respectively The couplings hjj, h D and h E are 3x3 Yukawa matrices, and 
\i and €i are parameters with units of mass. 

Supersymmetry breaking is parameterized by the standard set of soft supersymmetry 
breaking terms 

V soft = MfQTQ* + MfU*U,j + M'ifirhj + MfL«*L a J + MfR*Rj 
+m 2 H H a l *H a 1 +m 2 H2 H?H a 2 
- [±M 3 A 3 A 3 + \M\ 2 \ 2 + |M'AiAi + h.c] 



+£ab 



-BuHZHl + B&LSH^ , (1) 

Note that, in the presence of soft supersymmetry breaking terms the bilinear terms 
proportional to the e« can not be rotated away except for the very special case Bi = B 
and only if the scalar masses are adjusted in a special way. Such 'fine-tuned' assumptions 
at a low scale are, from our point of view, very unnatural. If realized at a high scale such 
as the unification or GUT-scale, then the trilinear R-parity breaking couplings introduced 
as a result of the rotation would re-introduce bilinear terms at the electroweak scale due 
to the structure of the corresponding RGEs |3"9fl . In contrast, bilinear terms are closed 
under RGE evolution from the high scale to the electroweak scale |^D| . 



In order to compare the t\ — > br decay mode with the 3-body MSSM modes it is 
sufficient for us to consider a 1-generation ]fi p model. Note however, that for the detailed 
connection of stop decays with neutrino physics we must consider the complete model, 
and we will do so in a second step when we discuss the connection with the solar neutrino 
mixing. Finally, notice that we also allow for R-parity-conserving Flavour Changing Neu- 
tral Currents (FCNC) effects, such as the process t\ — > c xl involving the three generations 
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of quarks in both cases. 

Our 1-generation model is specified by the superpotential ET 



Til 



W = h t Q 3 U 3 H 2 + h b Q i D 3 H 1 + h T L 3 R 3 H x - fiHiH 2 + e 3 L 3 H 2 



(2) 



For simplicity, in the remaining part of this section we adopt the 1-generation model 
when presenting the formulas for the mass matrices which are needed in the subsequent 
sections. The complete formulas for the 3-generation case are given in the second paper of 
[ |l3|j and have been used whenever required. The electroweak symmetry is broken when the 
two Higgs doublets, Hi and H 2 , and left slepton doublet L 3 acquire vacuum expectation 
values. We introduce the notation: 

■j= 2 [e 1 + v 1 +i ( p 1 ] 



Hi 



, H 2 



■±= 2 [v? + v 3 + iv I T ) 



The mass of W is given by m 



2 

W 



\g 2 v 2 , where v 2 = v\ + v\ + v 2 ~ (246 GeV) 2 . We 



define tan/? = v 2 j ylv^ + v 2 ). In addition to the above MSSM parameters, our model 
contains three new parameters, e 3 , v 3 and -B3, of which only two are independent, because 



there is an additional tad-pole equation [31]. These may be chosen as e 3 and v 3 . 
The stop mass matrix is given by 



Mr 



M 2 + \v 2 h t 2 + A UL ^ {v 2 A t - fivi + e 3 v 3 ) 



^ (v 2 A t - fMVi + e 3 v 3 ) M 2 + \v 2 h t 2 + A UR 



1 „Vi 



with A UL = \{g 2 - \g' 2 ) (v 2 - v 2 + v 2 ) and A UR = \g 
for the sbottoms is given by 



v 2 ). The mass matrix 



M 6 2 



M Q + ¥\ h l + A DL ^(VlA D -IIV 2 ) 

^ (viA D - nv 2 ) M 2 D + \v 2 h 2 + A DR 

where A DL = -\{g 2 + \g' 2 ) (vj - v 2 + v 2 ), A DR = -^g' 2 (vj - v 2 + v 2 ). The mass 
eigenstates are obtained by (q = t,b): 







cos 6q sin 6g 




Ql 






h 




— sin 6q COS 6q 








Qr 
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with 



cos% 



-M| 
gi2 



V /(^ 11 --|) 2 + (M| 12 ) 2 ' 



sin 8g 



M L - K 



(M~ 2 -m2) 2 + (M~ 2 



2 \2 

qi2' 



The bilinear term in Eq. (Q) leads to a mixing between the charginos and the r-lepton 
which induces the decay t\ — > br. The mass matrix is given by 



Mr 



M 







j^gv 3 -e 3 ^Mi 



As in the MSSM, the chargino mass matrix is diagonalized by two rotation matrices U 
and V 



ITMcV- 1 



m-± 

















m xf 



The lightest eigenstate of this mass matrix must be the tau lepton {j ± = xf) and so the 



mass is constrained to be 1.77703^Q2g GeVfH 



In our model, the one of the three neutrinos acquires mass at the tree level due to 



a mixing between the neutralino sector and one of the neutrinos $T\\ , |12| , f46| . The neu- 
tralino / neutrino mass matrix is 



M 



N 



M' 


-\g'v! 

\g'v 2 

-\g'vz 



o 

M 

\gvi 

-\gvi 

\g^ 



\g'vi \g'v 2 -\g'v z 



\gv x 



ft 







-\gv 2 



o 

£3 



1^3 



£3 





and M' is the U(l) gaugino soft mass. This neutralino/neutrino mass matrix is diagonal- 
ized by a 5 x 5 rotation matrix N' such that 



N'*M N N 



&&g{mx 1 ' m & m xp m x° i ' m x J 
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with m V3 = m^o. Note that in |26|, a different basis for the neutralinos was used. 



Assuming small R-parity violating couplings the tree-level neutrino mass is approximately 
given by 

(g 2 M' + g ,2 M)ii' 2 



m V3 « 



AMM'fi' 2 - 2(g 2 M' + g' 2 M)fi'v 2 v' 1 cos£ 



-v 1 sm 4 



(3) 



with 



sin£ 



e 3 v x + /iv 3 



/i 2 + e 



3 ! 



v I + v I 



Notice that in mSUGRA models with bilinear R-parity violation [32| m U3 is calculable 
through the RGE evolution and one finds in this case cancellations up to two orders 
of magnitude for the combination A 3 = 63^1 + /!t> 3 . In general models the smallness of 
m U3 requires relatively small e 3 values from the start as might arise, for example, in the 



models considered in ref. [21|. The remaining two neutrinos acquire mass radiatively. 
Rigorous quantitative results were given in the second paper in ref. | l3f . Typically they 
are hierarchically lighter than the heaviest neutrino, whose mass arises at the tree-level. 
This way one accounts for the observed hierarchy between the solar and the atmospheric 
neutrino mass scales. 

Similarly, the Higgs bosons mix with charged sleptons and the real (imaginary) parts 
of the sneutrino mix with the scalar (pseudoscalar) Higgs bosons. We denote the scalar 
bosons by the pseudo-scalar by Pf and the charged bosons by Sf. The relevant 
formulas have all been given e.g. in [^, |57[| . However, since we will take one of the 
pseudoscalar mass eigenvalues as input, it is worth briefly repeating here the discussion 
of the pseudoscalar bosons masses. The pseudo-scalar mass matrix is given by: 



-m ~Bse 3 /ie 3 ^ - B 3 e 3 f^ 



t':i 



(4) 



As expected, this matrix has zero determinant, since the neutral Goldstone boson eaten 
by the Z is one of the corresponding states. Therefore, the masses of the two physical 



states are given by the formula: 

m 2 , 3 = ^TrM ± (TrM) 2 - 4(M n M 22 - M 2 2 + M n M 33 - M 2 3 + M 22 M 33 - M| 3 ). 

(5) 

Therefore we can easily take one these masses as input and calculate B/i from it using 

-m 4 p0 v 1 v 2 v 3 + B 3 e 3 fiv 3 (v 2 + v\ + i>£) + e 3 m 2 p0 [nv 2 (vf + vf)) - B 3 v x {vl + vf) 



-m 2 p o{vf + vf)v 3 + e 3 (^t)i - B 3 v 2 )(vf + v\ + v\ 



B 3 is obtained from the minimum equation for given e 3 and v 3 [[35 1. 



III. NUMERICAL ANALYSES 

In this section we present our numerical results for the branching ratios of the lighter 
stop t\. Here we consider scenarios where all two-body decays induced at tree-level are 
kinematically forbidden except the b l + decays. Before going into detail it is useful to have 
some approximate formulas at hand 



w ~ j \ g 2 \U 32 \ 2 h£ cosj m h g 2 \e 3 \ 2 h 2 b cos j m h 

Tih^br)^ * « f (6) 

107T l07T|/i| Z 

/ m 2 Q \ 2 

m^cx^^Fht^coset-dAsmdtfflm^ I 1 - -Ji J , (7) 

where F = £ (log(m Gl/T /m z )if cb ^/16vr 2 ) 2 ~ 6 x 10" 7 , f L = v/2(tan 6 W N n + 3iV 12 )/6 
and the parameter <5 m 2 is given by 

_ M| + M|) + + A 2 , 

mi o 2 

m~ — m~ 

For the minimal SUGRA models one finds 5 m 2 = 0(1) which is basically independent 
of the initial conditions due to the mo dependence both in the numerator and in the 
denominator. Finally we have 

m t (A b + \A t ) 

d A — o 2 

m~ — VTir 
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The complete formulas are given in |!T], [32]], while for the three-body decays they are 
given in Appendix |A]. They reduce to the ones given in |27j for vanishing R-parity 
violation parameters. 

We have fixed the parameters as in |27] to avoid colour breaking minima, while in 
the top squark sector we have used mf , cos6^, tan/3, and fi as input parameters. For 
the sbottom sector we have fixed Mq , and Ab as input parameters whereas for the 



charged scalars we took m P o, M^, M^, and A T as input [|5T]. In addition we have chosen 
the R-parity violating parameters €3 and v 3 in such a way that the heaviest neutrino mass 
is fixed with the help of Eq. ([3]). For simplicity, we have also assumed that the soft SUSY 
breaking parameters are equal for all generations. 

In order to get a feeling for minimal size of branching ratios that can be measured let 
us first shortly discuss the expected size for the direct production of light stops at future 
colliders. One expects for example at the LHC a production cross section of ~ 35 pb 
for 220 GeV stop mass. Therefore, once the full luminosity has been reached, one has to 
expect approximately 3.5 10 6 events per year. The corresponding stop production cross 
section at a future e + e~ linear collider of 800 c.m.s. energy is of O(10 — lOOfb) ||23|| . For 
an integrated luminosity of 500 fb per year one can expect O(10 4 ) stop pairs per year. 
This implies that branching ratio as low as 10~ 3 can in principle be measured. 

We consider first the simplest case of one generation model which, as already men- 
tioned, is sufficient to describe the relative importance of the t\ — > t + b decay mode 
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relative to the possible 3-body decay modes 

t~i - W+bx\ 

ti - s+bxi 

ti -> S£bv 3 

ti^bSfr + , (8) 
ti -> &P°r + 
ti -> 6Zf i/j, 

ti ^ bDi l + {I — e, fj) . 

In general the important final states are those that conserve R-parity. For example, for 
the case of decays involving S£ the most important are those in which the scalars mainly 
a stau. Due to the fact that the existing bounds on the MSSM sneutrinos are below 
100 GeV there exists the possibility that the sneutrino has nearly the same mass as one 
of the Higgs boson. Similarly it could be that the charged MSSM boson has nearly the 
same mass as one of the staus. This implies large mixing effects even for small R-parity 



breaking parameters |35], fI7|| . We therefore have used the complete formulas for the 
3-body decay modes which are presented in the Appendix. The latter include R-parity 
violating decays such as t\ —>■ W + bi > 3 . In addition to the above mentioned decays there 
is also ti —>■ b Z° t + . This decay mode is kinematically suppressed compared to t\ — > br + 
and there is no possible enhancement due to a mixing with an R-parity conserving final 
state. Therefore it can be safely neglected. 

In Fig. |l] we show the branching ratios for the ti as a function of cos 6$ in different 
scenarios. In order to calculate the partial width for the decay ti — > cx\ we have taken 



the formula given in ref. p5[ . According to the analysis performed in |32|| , where the 
full calculation was done in the mSUGRA scenario, the result obtained with the present 
approximation should be taken as one upper bound. This implies also that the shown 
branching ratio for t\ — > br + can be viewed as a lower bound. The parameters and 
physical quantities used in Fig. [1] are given in Tab. Q. For the case of Fig. |l](a) we have 
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fixed in addition the R-parity violating parameters such that m v ~ 1 eV. With this choice 
of parameters S®, S®, P° , and S± are mainly the MSSM Higgs-bosons whereas S®, P%, 
5*2", and 5*3" are mainly the MSSM sleptons of the third generation. In the plot we show 
the various branching ratios of the lighter stop summing up those branching ratios for the 
decays into sleptons that give the same final state, for example: 

ti -> b u e e\ -»• b e + v e \\ , t 1 ^be + v e -> b e + v e \\ . 

The branching ratios for decays into jiL or are practically the same as those into 
cl or v e . Note that the energy spectrum of the leptons in the final will be somewhat 
different depending on whether the scalar in the intermediate step is charged or neutral. 
This offers in principle the possibility of determining the branching ratios of the different 
decay chains even if the final state topology is common. Note, that states containing 
scalars or neutralinos will lead to additional jet and/or lepton multiplicities absent in the 
MSSM. 

In Fig. 0(b) the slepton mass parameters are chosen such that decays into scalars 
are kinematically forbidden. Here we display the channels t\ — > bW + Xi, ^1 ~~ * 
and t\ — > cxi- The remaining modes, such as t\ — > 6S , °r + , turn out to be completely 
negligible. In both cases, with and without sleptons in the final state, one can see that 
in general the three body mode t\ bW + Xi, dominates except for a somewhat narrow 
range of negative cosd^. However, the branching ratio for t\ — ► br + is above 0.1% for 
most values of | cos#j| implying the observability of this mode. Most importantly, note 
that even in the parameter ranges where the three-body decay mode is dominant, its 
resulting signature is rather different from that of the MSSM due to the fact the lightest 
neutralino decays into SM-fermions, leading to enhanced jet and/or lepton multiplicities, 



as discussed in detail in |37|, |38|]. In the remaining part of this section we assume that 
3-body decays into scalars are kinematically forbidden. 

In Fig. U(c) the R-parity violating parameters are fixed in such a way that the heaviest 
neutrino mass is in the range suggested by the oscillation interpretation of the atmospheric 



neutrino anomaly [14 
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FIG. 1: Branching ratios for the t\ as a function of cos% for different scenarios. We have 
fixed in a) m Vz = leV, b) m V3 = leV, > 225 GeV, M L > 225 GeV, c) m Vz = 0.06 eV, 
Mg > 225 GeV, > 225 GeV, d) Branching ratios for the t\ as a function of cos 0% for 
tan/3 = 3. m U3 = 0.06 eV, > 225 GeV, M L > 225 GeV. All the other inputs are given in 
Table |. 
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Input: 


tan P - 


= 6 


/" = 


500 GeV 


M = 


250 GeV 




D 


370 GeV 


Ma 
Q 


= 340 GeV 


A b = 


= 150 GeV 




Mp = 


210 GeV 


M f 


= 210 GeV 


A T = 


= 150 GeV 






220 GeV 


cos 6 


t = ~ - 8 


Tfl pO 


= 300 GeV 


Calculated 


m-o = 


122 GeV 


m~+ 

V J 
A 1 


= 234 GeV 


m~+ 


= 519 GeV 




777.1 = 


334 GeV 


"% 2 


= 381 GeV 


cos 9 


I = 0.879 




m s o = 


107 GeV 


rriao 


= 200 GeV 


m s o 


= 302 GeV 




m p o = 

^2 


= 200 GeV 


^3 


= 300 GeV 








m Q - = 


= 203 GeV 


Trio- 
rs 


= 226 GeV 


rrio- 


= 311 GeV 




mg L = 


215 GeV 




= m 9ii = 200 GeV 







TABLE I: Input parameters and resulting quantities used in Fig. |l]. 

In Fig. |I|(d) we show the same scenario as in Fig. |I|(c) but for tan f3 = 3. The branching 
ratio into br now increases, whereas the branching ratio into cxi decreases. This is easily 
understood by inspecting Eqs. (j^) and (0). Indeed for the br case the partial width is 
proportional to h%, whereas for cxi it is proportional to h\. This implies that the partial 
width for t\ — > cx\ grows faster with tan/5 than the width for t\ — > br. This is also 
demonstrated in Fig. ^ where we show the tan /3 dependence of the branching ratio for 
the decay of t\ into br + for several values of the neutrino mass. For m U3 = 0.06 eV the 
B(t\ — > br) is still above 0.1% if tan/3 is not too large, as favored by the explanation 
of the neutrino anomalies in this model fl3| . As seen from the figure, the the t\ — ► br + 
branching ratio is also somewhat correlated to the U3 mass. Should one add a sterile 



neutrino to the model f48| , then the neutrino state z/3 could in principle be heavier than 
assumed above, favoring t\ — > r + b decay mode. 

Let us now turn to the general three neutrinos case. There are new features that arise 
in this opposed to the 1-generation case considered so far. In this model the 

solution to the present neutrino anomalies implies that all the £j are of the same order of 
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B( fj -> b x ) 



m .=lOO eV 




FIG. 2: Branching ratios for t\ decays for m fl = 220 GeV, /j = 500 GeV, M = 240 GeV, and 
m u = 100, 1 and 0.06 eV. The branching ratios are shown as a function of tan (3. (cos 9f = —0.8) 



magnitude [13[ 



Two furhter important results of |T3[ are that the atmospheric neutrino angle is con- 
trolled by the ratio fav d+^v 2) / (e^v d+H-v 3) and that the solar mixing angle is controlled by 
( e i/ e 2) 2 - One can get approximate formulas for the decay widths t\ — > be + and t\ — > b fi + 
similar to Eq. ([]) by replacing €3 by ei^. This implies that (i) The decays into be + and 
6 /i + are as important as the decay into br + . (ii) The decays t\ — > 6e + and ti b fi + are 
related with the solar mixing angle. Moreover, we find that X]z= e u r ^(^1 ~~ * ^ + ) * n ^ e 
3-generation model is nearly equal to T(ti — > br + ) in the 1-generation model provided 
that J^i=i e ? i s identified to e 2 in the 1-generation model. 

In Fig. |H we show the ratio of B(ii — > 6e + )/B(t! — > versus (ei/e2) 2 for different 
values of cos 6~ t . For definiteness we have fixed the heaviest neutrino mass at the best-fit 
value indicated by the atmospheric neutrino anomaly. One can see that the dependence is 
nearly linear even for rather small cos6^. For | cos^| < 10~ 2 the approximation in Eq. (0) 



breaks down and additional pieces dependent on sin^ |3T| , p^] become important, leading 
to the non-linear dependence. One sees from the figure that, as long as cos# t - > 1CT 2 
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B(fi -» be+)/B(ti ->■ 
1 

0.1 

10" 2 



10" 



10" 









> 0.1 // 




y/\ 2 
// 10~ 2 


: io~ 3 / 








0.1 









10" 4 10" 3 10" 2 0.1 1 

FIG. 3: Ratio of branching ratios: B(ii — ► 6e + )/B(ti — ► as a function of (ei/e2) 2 f° r 

m h = 220 GeV, fj, = 500 GeV, M = 240 GeV; | cos%| > 0.1, 0.01, 10~ 3 , m Vz = 0.6 eV. 

there is a good degree of correlation between the branching ratios into B(t\ —>■ be + ) and 
B(t\ — > and the ratio {t\/t2) 2 - Thus by measuring these branchings one will get 

information on the solar neutrino mixing, since tan 2 8 so i is proportional to (ei/e2) 2 |13 



which makes it a rather important quantity. For the so-called small mixing angle or SMA 
solution of the solar neutrino problem we expect t\ —>■ e + b to be negligible. In contrast, for 



the large mixing angle type solutions (LMA, LOW and QVAC, see ref. JTJ] and references 
therein) we expect all ti — Z + 6 decays to have comparable rates. As a result in this 
model one can directly test the solution to the solar neutrino problem against the lighter 
stop decay pattern. This is also complementary to the case of neutralino decays considered 
in ||38|| . In that case the sensitivity is mainly to atmospheric mixing, as opposed to solar 
mixing. Testing the latter in neutralino decays at a collider experiment requires more 
detailed information on the complete spectrum to test the solar angle In contrast 
we have obtained here a rather neat connection of stop decays with the solar neutrino 
physics. 

Note, that this result is much more general than the scenarios discussed in this paper. 
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It is of particular importance in scenarios where only the R-parity violating decays and 
the decay into Xi c are present |TT], Similarly, the other ratios of the final states bl + 
are proportional to the square of the ratio of corresponding 6j provided that cos Q~ t is not 
too small. 



IV. CONCLUSIONS 



We have studied the phenomenology of the lightest stop in scenarios where R-parity 
violating decays such as t\ — > b r + compete with three-body decays. We have found that 
for < 250 GeV there are regions of parameter where t\ — > b r + is an important decay 
mode if not the most important one. This implies that there exists the possibility of full 
stop mass reconstruction from r + t~ bb final states, favoring the prospects for its discov- 
ery. In contrast, in the MSSM the discovery of the lightest stop might not be possible at 
the LHC within this mass range. This implies that it is important to take into account 
this new decay mode when designing the stop search strategies at a future e + e~ Linear 
Collider. Spontaneously and bilinearly broken R-parity violation also imply additional 
leptons and/or jets in stop cascade decays. Looking at the three generation model the 
decays into t\ — » bl + imply the possibility of probing ef/ e| and thus the solar mixing 
angle. This complements information which can be obtained using neutralino decays. In 
the latter case the sensitivity is mainly to the atmospheric mixing, as opposed to solar 
mixing. In this model neutralino decays is ideal to test the atmospheric anomaly at a col- 
lider experiment, while stop decays provide neat complementary information on the solar 
mixing angle. Obtaining solar mixing information from neutralino decays would require 
more detailed knowledge on the supersymmetric spectrum, since it would be involved in 
the relevant loop calculations of the solar neutrino mass scale and mixing angle. By com- 
bining the two one can probe the parameters associated with both solar and atmospheric 
neutrino anomalies at collider experiments. 
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APPENDIX 

In this set of appendixes we present the formulas for the lighter stop decay widths and 
couplings used in the paper and omitted in previous sections. 

APPENDIX A: FORMULAS FOR THE THREE-BODY DECAY WIDTHS 
1. The width r(ti -► W + bXi) 



16 7rm? sin 4 ^ 



w 




+ 



+ Gl t + Gl i + Gf t+ Gl + Gl) 



(m b +m~o) 2 



i 
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The integrals Jtttst are: 
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with % — 0, 1, 2. Their integration range is given by 
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where s = (p^ — £>vk) 2 and £ = (p^ — p t ) 2 are the usual Mandelstam variables. Note, that 
— T~+m~+ appears in the entries of the integrals G^+i. and G^ +f because the chargino is 



exchanged in the w-channel in our convention. The coefficients are given by (no sum upon 
repeated index): 
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G tt = (d ikl + d ik2 s)J°(m 2 , T t m t ) + (d ik3 + d ik4 s)J^(m 2 , T t m t ) + d ik4 J 2 (m 2 , T t m t ) , 



z 

G tb = ^[(e ikll + e ikl2 s)J° t (ml i ,r i m ii ,m 2 ,r t m t )+ e ikl3 J^(m? , T~ b m~ bi , m 2 , V t m t ) 
i=i 



Gz 



^\(s, m|, m 2 s +)X(s, m|o, m 2 ) 



x 



+ Re 



(/ifc31 + fik^s) 



{s - m 2 bi + ir bi m bi )(s - m 2 ^ - iT b m b2 ) 



^ is — m~ ) 2 + r~ m~ 

l=i K h> h b t 

Their integration range is given by 

+ ml + m 2 s+ + m|o - s (to? - to| + )(to| - to|) 



m 



ti 



t max 
min 



2s 



± 



^\(s,m 2 i ,m 2 s +)\(s,m 2 c o,m 2 ) 



2s 



where s = (p^ — Ps+) 2 an d t = (p^ —pt) 2 are the usual Mandelstam variables. Note, that 
— T-+m-+ appears in the entries of the integrals G^ +b . and G^+t because the chargino is 
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exchanged in the w-channel in our convention. The coefficients are given by (no sum upon 
repeated index): 



-4 fcyijQ^ Q^Wri^o (ml + m 2 .+ + m|o + m\ + m 2 ^ 
-2 Qt jk 'Q* k [(k^) 2 + (/*,.) 2 ) m ftj+ (2 m 2 + m|o + m 2 +) 
"2 fc{,-4- ((^') 2 + (Q* k T) m b m n (ml + 2 m|o + m|) 



- (^) 2 (Q^) 2 + (C) 2 (^V) 2 



2 1 2 
m 6 + m -0 



: ( m 2 + m 2 + 



2 1 2 



((*L-) 2 (OSV) 2 + (il-) 2 (QS*) 2 ) ™~+ H + 4?) , 



+2 fc^L- m 3k 'Y + (Q? jk 7) m b m~+ + 2 Q^Q*' ((^) 2 + (^) 2 ) 

+ ((*L-) 2 (oSV) 2 + (iloW) 2 ) ™ 2 + , 



A.j A. -i 



+ ((*L-) 2 (QS*') 2 + (40 2 (QJV) 2 ) (2^ + 2m| ? + ml, + ml) 
+2 k\^ ((QJV) 2 + (Qg*') 2 ) ™^x + » 
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&ijk4 — 
&i4fcl = 



(Qijk ) (hj) (Qijk ) ) 
-4 (^AQf^UV + k\ 2 l[ 1 Qiik'Q?2k) ™ b mto (m|, + m|+ + m 2 + m|) 



-4 (M2*iiQ,VQ2fc' + KA2Qiik'Qf2k) m * m n m xt m xi 

—2 ( llil^Qiik Qi2k + kiikuQak Qf 



'%2k 



X 



(ml + m|o) + (m|o + m?) (m 2 + m 2 + ) 



2 ( k^l^Qiik Q^k + k^nQak Qa,k ) m bm- 



Xi 



tl 



2ml 



-2 fe^V^V + ?iAQwfQm) ™*+™x+ + m ' 



+4 ^n4 2 Qi2fc Qflk + ^L^llQflfc Qilfc ) m 6 m xJ 

+2 (fc'nfc^'Ofi*' + ™x°™*+ 

+2 {l^A^Qiik Qa.k + ^11^12 Qufc'QiL') ( m &+ m |°) 

+2 ^fciiii2 Qafc Q^fc + ^L^iiQufc QiL ) va^m^r 

+2 {^AiQiik Q%,k + Mi^Qufc'QiL') m ^> m xi 

+2 (klikl 2 Qf lk 'Qi2k + ^L^afc'QiL') m xt m xi ' 
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2 (jn^nQiak Qak + Ki^uQiik Q 



i2k 



+4 
+2 
+2 
+2 
+2 



m^om~+ 



^11^12Qi2k Qak + k^llQilk Qi2k 



k u k 12 Q i2 k Qak + hih2Qak Qv2k 



it it f^L if^L I _| it jLt <oR / 

jLt 7* f\L I f}L I | Kt 7t /^i? / 

ft ll ( '12Vilfc Vj2fc "r /t 12 t llVilfc Vj2fc 

7t it f)L lr\L I , Kt 7t ,o,R /^i? / 
,i 12 t llVilfc Vj2fe ^ ,i ll t 12VjlA: Vj2fe 



2 + 2 m|o + m|+ + m| 



m b m~+ 



The coefficients a^z are obtained from a^/ by replacing: l\ 2 — > Z* 3 /c* 2 — > /c* 3 

Qm Qak ~^ Qi3k' m x+ ~^ m x+ an d the coefficients a^ki are obtained from a^ki 

replacing: a i4kl -> a i6fe/ : Z' n -> Zj 3 fc' n -> k{ 3 Qf lk ' -> Qf 3fe ' Q* fc ' -> Q^' m-+ 



m x 3 + 



v 



bukljQfjk'm^+mbmt 



R 



5± 



«1 



2 2 2 2 "\ ~^k2 2/2 2 r\ 2 "\ 

m s+ m ti ~~ m b m x° J — m b yT^b + m t 1 + 2 m x 9 J 



«1 



X 



2 / 2 , 2 2 

m 6 ( m i<7 + + m fi — m-q — m 



■5± 



,2 f__2 , .__2 „2 „„2\ R-kl 



2 / o 2 1 2 1 2 \ 

ra t 12 m fe + m -p + m 5 + I 



rL k2 
V 2 



rn 



d5± p5± 
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+b\ i k\ j QfJm^m b m t 



x 



2,2 2 

m fi + m s + —m b - 



—o\ik\jQ^ k rrib 



x 



m l m x ~ m \ m \+ 



R 



t>2 



^1 



5± 
fcl 



(m|o-m|)^+m t 2 (m2 + m|o) 



fc2 



V2 



5± 

fcl 



pS* dS± 

-ft 1.0 , ? , + ^ f? < -ft 



^2 ^2 



2^ 



Hjk3 



b\ i k t lj Qf jk 'm~+m b m t 



R 



fc2 



5 

+«H A; L' ( 5 J L 1 'fc ,m y+ m fe m >?? 



v 2 J J 1 f 2 



5± 
fc2 



+ m|+ + m| + m 2 , 



fc2 



5± 



^2 



2 rl fcl 



«1 



fcl 



blikljQfjk'm^ombmt 



"u^ij^ijk •• v xJ" v °" v xt ,.. 



5± 
fcl 



fc2 



5± 
fcl 



+a t li k\jQ^ k rrib 



2 i 2 i 2 i 2 \ 

m fe + m-o + m s + + m f l 



R s 

2 \ ""-fcl 



2 rl fc2 



^1 



^2 



2^ 



blAjQijkmt 



R k2 

V2 



+ O-likljQ^jk m b 



R kl 

V2 
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Cijkll 



2C k 



(bllljQijk + m b (m 2 h + m\ + 2m|o) 



+ 
+ 



Cijkl2 — 



-2C- t \ 



blAjQijk + a U^ljQ^k ) m b + [ a ulljQijk + ^li^ljQfjk 



+ \biik\jQijk + a iil\jQijk 



m~ 



x 



m^o 



AjklZ 



{tf,iAjQijk 



+ 




Hk2 



2m 2 t 

IF 



"li, 



R 



A; 2 



V2 



+ m 



R ki 

Vi 
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dm = -JU ((4) 2 + (b\ t ) 2 ) 



m\m 2 



2a u b H m t m^ 



mil 2 + 



R 



fei 



R s± 
n ki 

Vl 



Vl 
2 



R s± 
- Ill 2 1 "'■ - 

' ' ^2 



V 2 



ik4 



R 



5± 
fcl 



t ^22 I rL fe2 



«1 



^2 



2^ 



fa-*) 



R 



A; 2 



+b\$ u m h m t 



V2 



2m 



2- ,L fcl 



«1 



H.-^)^-( 



5± 



fc2 + mg 



2 / 2 , 2 \ -^fc2 , 2 { 2 2 \ 

m t (m b + m^ — +mb [m-o - m f J 



5± 
fcl 



-6fci a i* m 6 m v? 



2m 2 t ?S-+(m 2 b - ml + ) ^ 



5± 



_ '-V^ >/ ( n b n i m 2 R h2 ,J> U "AJ 

e^2 — C figj ( a H a H m t 



s± - - Rt 
+ b b ki b\-m h m t - 1 



v 2 



Zikl3 



2V2 



R 



A; 2 



(.6 !,* 



5± 
fc2 



^2 



V2 



R 



5± 



- - R s± - - 
+ 4^™%-^ + b h ki a\^n h m tl ^ 



fikll — 



:4) 2 + 



m 2 b + m|o) +4a^m^o] , 



/*» = (c^) 2 ((4) 2 + (6Si) 2 ) > 
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ik31 



-2 C-\ C~l 
ti&i tib2 



4<4 + b b u b b 2i ) [mi + m\, J + 2 f a° u b° 2i + b° u a° 2i J m 6 m^o 



,6 ub 



/*» = 2C$C$ 2 {a^A + b\^ 



3. The width r(*i -»• S% b xj 



The decay width is given by 



nti^sibxt 



a 



16 7rm? sin 4 6-w 
ti 



Kr m so) 2 



(m b +m +) 2 
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with 

3 

G f+x+ = [( a ^ fcl + a ijk2s)J b 3 (m 2 ii + m|o + ml + mi+ - ml+ - s, F-+m~+) 

3=1 

+(dijk3 + aiAjkS)Jl{m 2 7 + m| + m 2 b + m~+ - - s, r~+m~+) 

+ a-ijki Jb( m i + m s° + m 2 b + m l+ - m l+ ~ s ? r y+ m y+) 
+(a mi + a m2 s) J° (< + m|o + + + - m|+ - s , r~+m-+ 

> m| + m|o + m 2 b + mj+ - + - s , r~+TO~+) 
+ (a^ 3 + a^s)Jl b {ml + m|o + m 2 b + m|+ - m| f - s, r~+m,+ 

> m| + m|o + + + - + - s , T-+m- 2+ ) 
+ a iAkA Jl h {m 2 h + m|o + m 2 b + - m;?+ - s, T^+m^+ 

, m? + m 2 s o +m 2 b + + - mj+ - s, r^ + m-+) 
+ (ai 5 fei + ai<ok2s)J bb {m 2 h + m|o + + - - s, T-+m-+ 

. m l + m l° + m l + m \+ - m \+ - s > V xi m xt) 

+ (a i5k3 + a ibM s)Jl b (m 2 h + m|o + m 2 b + - - s, T^+m-+ 

> m !i + m l° + m 6 + m l+ - m l+ - s > T xi m xt"> 

+ a i5k4 Jl h {m\ + m|o + + m|+ - - s, T-+m-+ 
> m| + m|o + + + - m | + - s , r-+m~ 3+ ) 



+ (aj 6 fei + am2s)J bb (m 2 i + m 2 s o + m 2 b 


+ m x+ 


— m 2 + 

xt 


-^ T xt m xi 


2 , 2 , 2 

, m h + m s o + m b 


+ ™x+ 


— m 2 + 

xt 


- s ^xt m xt^ 


+ (a i6 k3 + aiQ M s)Jl b (m 2 h + m|o + m 2 b 


+ ™x+ 


— m 2 + 

xt 


- s > T xi m xt 


2 i 2 i 2 

, m^ + m s o + m b 


+ m 2 + 

xt 


— m 2 + 

xt 


-^ T xt m xt) 


+ a m4 Jl h {m 2 h + m 2 s o + m 2 b + ml+ - 


m 2 + — 

xi 




, m\ + m|o + + + - 


m 2 + — 

xi 
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r,s° 



^2 ( bi i kl + b ijk2-s)Jb b (m 2 ii + m|o + m 2 b + mi + - mi + - s, -T^+m^+,m 2 b , Y b m b ) 



+(b ijk s + b ijki s)Jl h {ml + m|o + m 2 + m\+ - m\+ - s, -T^m-pm 2 , T b m b ) 
+ b ijk A J bb {m 2 + m 2 s o +m 2 b + m 2 + - m 2 + - s, -r~+m~+,m 2 b , T b m b ) , 

1 k A,i Aj J 



G x+t 



z 

^ |^(cijfcn + Cijkns) 
1=1 

* J st( m h T t i m i[ ,m 2 i + m 2 s0 +m 2 b + m 2 + - m 2 + - s, -r~+m~+) 

ft A. j J 3 J 

+ c ijk i3 Jlt( m l > r*," 1 *, > m l + m h + m l + m l+ - m l+ - s > - r v +m x + ) ' 

ft ^-i •Aj J J J 



Gf t ° = (djfci + d ik2 s)J b (m 2 b , Y b m b ) + (d ife3 + d iM s)J b (m 2 b , Y b m b ) + d ifc4 Jjf (mjj, r fe m b ) , 



70/2 



x 



/ A(s, m\ , m| )A(s, m 2 + ,m 2 b ) 



(fikll + fikl2S) 



i=i K ti> ti ti 
The integrals Jt'tt,st are: 



+ Re 



(/ifc31 + fik32$) 



s - m? + ir t - l m t - 1 )(s - m| - ir t ~ 2 m t ~ 2 ) 



t 2 tJ 't2"^2/ 



tmax 

dt 



(t - m 2 f + mfr? ' 



J b l b (m l ,m l T 1 ,m 2 ,m 2 T 2 ) = Re /rft 



{t — m\ + imiFi)(t — m 2 — ?m 2 r 2 ) 



J l st (m 1: mi ri,m 2 , ra 2 r 2 ) = Re 



dt- 



s — + im\Yi J [t — m\ — im 2 Y 2 ) 
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with % — 0, 1, 2. Their integration range is given by 

to? + m 2 + to| + to 2 + — s (m~ — to| )(to 2 + — to 2 ) 



t 



max 
min 



2s 

/ A(s, to? , mj„) A(s, TO 2 , TO 2 ) 

2l ' 

where s = (p^ — p 5 o) 2 and t = (p^ —pt) 2 are the usual Mandelstam variables. Note, that 
— T-+TO-+ appears in the entries of the integrals and G?+ t because the chargino is 

exchanged in the -u-channel in our convention. The coefficients are given by (no sum upon 
repeated index): 



njfcl 



-4 k t lj l t lj lf- k kf- k m b m~+ (m 2 b + to 2 + + to 2 + + m? + m|o J 
-2 C*C {{k[jf + (/^.) 2 ) m- t m- t (2 to 2 + to 2 + + to|o) 

-2 fcyl- ((C) 2 + ( fc S) 2 ) H + 2 m *i + <) 
-((^•) 2 (0 2 +(4) 2 (C) 2 ) 

x (m 2 b + to 2 + j + (m 2 b + to|o) (to 2 + + to?) 
- ((^) 2 (© 2 + (4) 2 (C) 2 ) - 2 + H + ' 



a ijk2 = Ak[fygkgm b m^ + ((^) 2 (^° fc ) 2 + ( m l + m l+) 

+2 fey L- ((C) 2 + (C) 2 ) ™»™*+ + 2 CC ((*l-) 2 + (40 



\2 



m x+ m x+ 



+ ((feL) 2 (C) 2 + (^) 2 (^S) 2 H 2 ^ 



a^ 3 = 4feyyg° fc feg° fc TO 6 TO-+ + 2Zg° fe feg° ((fe^.) 2 + (Zy 2 ) m^+m^ 

+ ((^) 2 (C) 2 + (4) 2 (C) 2 ) ( 2 - 2 + <o + ml) 

+2 MA ((C) 2 + (C) 2 ) > 



a jjfc4 — (felj) 2 (Xgfc) 2 (^lj) 2 (^'fc) 2 ? 
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«> ia-j = -2 ^n^ 2 ^2°feS + k^kfiA) "i^ m xf ( m x+ + m ^ + 2m ') 

-4 + fc^n^iX) ™&™x+ ( m *+ + m ^ + m l + m l) 

s ) m b m x +r 

—2 ( k\ 1 k[ 2 kf 2 k^iik + ^n^i2^afc^fc ) m xt m xt 



x- 

- 1 ( k&li^A + m b m^m^m^ 



it it uS°iS° 



m/i + 4 m 2 s o +2 ml 



2 f /* /* A- 5 ° A- 5 ° 4- fr* k l 1 S ° l s ° } 

Z I l\\l\2H\kH2k + K ll hj 12 l ilk i i2k J 

(m 2 b 4 ml+^j 4 (m 2 .+ 4 ra|) (mg + m|o) 



-2 + fcL4€X) ™&™*+ H + m ^ + 2m lt) 

— 2 ^ii^i2^lfe^fc + ^11^12^ilk^i2kj m xX m xt ( m ^+ + m b 



S° ,5° 



2m 2 + 



&i4k2 — 2 ( ^n^i2^fc^afe + k\ik\ 2 k^ lk l^2k ) m x +m x 



„s° ,s° 



+4 
+2 
+2 
+2 
+2 
+2 



. t it > S° iS° , . t it .SO .S" \ 
^ll L l2^i2k L ilk + K l2 i ll K ilk i i2k J " lb " l xt 

h t h t j,S°jS° , it ft V.S°iS°\ 
^\1^12^i2k L ilk + L \\ L \2 K i\k L i2k J " L xt X.t 

ft f k S° k S° , 7 t > t iS° iS° \ ^ 2 I 
L \\ L \2 K i\k K i2k + K \\ K \2 L i\k l i2k ) \ rn b + " l 



Lt it . S° , , t it ,S° ,S° \ 

Rj \\ i \2 K i\k K i2k + K \2 i \\ l i\k l i2k J '''6''^+ 

li It i.S° r,S° , Kt It ,S° lS° \ 
^\2 L \\ K i\k K i2k + K ll i \2 i ilk i i2k J " lb ' n xt 

a a k s° k s° , it it is° is° \ 

K U K 12 K ilk K i2k + l ll l 12 l ilk l i2k J m xf X.t 



liikz — 2 (^{juk^Jiik + ^n^L^fc^fc) m xt m xt 
+4 + M^ii^ufc^fc) mbm xt 

-|-Z I ^\\^\2^i2k L i\k + i ll i 12^ilk i i2k J " l xt X.t 

+2 + M1M2SS) (2^5 + 2m^ + m|o + m|) 

+2 ^ii^i2^ifc^fc + ^L^ii^ifc^fe) mbm xt 
+2 r *■* /* A- 5 ° A- 5 ° 4- *■* /* / 5 ° / 5 ° ^ + 

I l^uhl^ilk^k + /t ll ( '12''ilfc''j2fc J rnbrn xt ' 
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Z I I>ui<i2^ilk nj i2k + Rj \\ K \2 L i\k L i2k 



The coefficients 0^5^ are obtained from by replacing: 7' 2 — > 7* 3 fc 



b 12 



A: 



13 



1.50 
ft i2fc 



m 



and the coefficients a^i are obtained from a^ki by replac- 



ing: a i4fc / — > a i6fc / : / 



11 



it ui 
l 13 ^ 



11 



k 



13 



ft ilfc ~ *■ ft i3fc l il7c 



l i3k 



m xt 



m xt 



b ijkl = 2l kliky^m^+ml (m\ - m 2 + ) 7f 6 ° - (m 2 b + m|+) k 



+k\ i l t i j lf- k m h (m z s oml - m£m*+ j 1% - m z b [m z b + + 2 m z .+ J 7c, 
4-7* 7* l S °m a 



7S 



"66 



2 7 2 1 2 

+k li k 1 jk i j k m^+m b 



m 



x 



m? I 7c 



) 



.5° 
"66 



2mf7 



2,5° 



66 



l bb - ( m |o + 2m 2 b + m 2 + ) k 



+k t u l\ j k^ k m b m~+m i 



(m 2 s o - mfj 



,2/ „S° 



1+ - m 2 ) 7cf 6 ° - m 2 (2 m 2 + m 2 + + m|o) 



66 



-l^kljkfjf.rrib 



x 



2 2 
m h m.+ 

° xj 



rn 2 _rii 



sg) 4° + ™6 (™6 + ™| + 2 m|+) / 



5° 
66 



-7* 7* k s m + 



ml+ - m| ) /cf h " + m 2 ( m 2 + m 2 + ) 



Xi 



bijk2 — 2-^ k^kljl^m^+mfk^ + k^lljlfj^.m^rribk^ 



J ii 7S 



2 7,5° 



,t 7* 75° , 



,5° 



+ 7* 7* / S 777 2 777 + 7 5 4- /V* /V* /V 5 %T7 4-TO 2 /V 5 " 
^ i li i lj i ijk" l b" l x+ i bb + ^li^lj^ijk" 1 ^ Ul bHb 



,2 7 S° 



t iJ uS°, 



,2 75° 
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ljhjk m b 



m. 



+ m|o + ml 
t iJ 



m 



2 W 



s° 



,2 U S° 



hi, m b^bb 



+l\ i k\jlfj k m^+mi } m^+k b s b k^k^kf^m^+m^k^ 



+kl i ll j kf° k m~ + m b m~ + l^ + AA^mim^k^ 



i it 1,5° 2„ 



,5° 



+l t li k\jkfj k rrib 



f 2 i 2 i 2 i 2 \ ; 5° 2 7 5° 

^m fe + m.+ + m 5 o + m t -J fc 66 - m h l 



bb 



%jki — 2 (kiilijlfjk m bliib + ^iikijkfj k mbk bb j , 



titi 



+ 



Cijkl2 



titi 



kjil\jlijk + luk[jkfj k ) nib + \ l\il\jlij k + k* H k[jkfj k j m-n 



( k\ik\jlfj k + ^Uljkijkj m x j 



Cijkl3 



liit jS° , ,t 7 t .5° A , / rt ,t ,5° , . t . t .5° \ 
K li L lj L ijk + Hi^lj^ijk J rn b + \ i li L lj L ijk + Hi K \j K ijk J " L x+ 



44 



di 



ikl 



(lifml (kSY + ikifml (if 



+ 



+2 ll^mbm^ 



2 ( 2 i 2 \ 

m 6 (^m fc + m .+ j 



m|o - m 2 ) ( m 2 ( /cf h ° ) + m 2 ( Zf h ° 



+2 ((4) 2 + (4) 2 )-^ (m 2 -m 2 + )j, 



2 mlml 



d 



ik2 



m\ 



hi, 



(k 



CO 



^3 = -|2((4) 2 + (A;y 2 )mX + 244^^ + 

fe) 2 + (4)X (if 



m. 



2 + (4°) 2 )+^(C) : 



m|o + m\ 



dikA — 



{llfml (4°) 2 + iklfml (IS] 



Zikll — 



-2^1^47^+ [(m|o 



m 2 ) lf-2m?£ 



6 ; i w> 



-/* /* 



2 / 2 , 2 \|S° 2 f 2 2 \ 1 S° 

m b [m b + m.+ J Z 66 + m 6 (m .+ - m fi ) k bb 



2m 2 /g' + (m 2 h -ml 



(m 2 b -m 2 s ^jk^ J, 



era, = -2 (4«C + 44^ 2 4°) , 



'tit; \ lukli m b m x+lbb + ^L^H m2 ^W) 



bb ) 
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ikll 



2 , 2 

m fe + m + 



+ AV li k t kl m h m- 



ftm = (^|) 2 ((4) 2 + (4J 2 ) • 



ifc31 



qO qO 

-2 C~ *~ cfi 

till 11*2 



/* /* _|_ U 1 h l 
l \i l 2i ^ h li h 2i 



m 



fik32 - 2C S i ^C S t ^ i2 [llfi 2i + k\ i k t 2i 



APPENDIX B: COUPLINGS 



Here we give the couplings that were used in sec. 
The (Ji-q'-xf couplings read then 



iz q o q 

in jni 



k 3 = nou 



with 



-v, 



htV n 



tV]2 



0\ 



-u jX 

h b U j2 



where Uij and Vij are the mixing matrices of the charginos [55|. 
For the couplings x^ ~ S°/P° — x ± we have 



L ijk 

k S ° 
^ijk 

P 

Hjk 

k P ° 
^ijk 



~ RkaVjiUis — -Rfc2^j'2^ii — R^iVjiUiz + h T (R% 3 Vj 3 U i2 — R kl Vj 3 U i3 ) 
- RkzUj-sVn - Rl 2 UjiV i2 - Rk l Uj 2 V il + h T {Rl 3 U j2 V i3 - R kl U j3 V i3 ) 
~ Rk3 VjiUi3 - Rk2 Vj2Un - Rfo VjiU i2 - h T (i?f 3 Vj 3 U i2 - Rfo Vj 3 U i3 ) 
~ Rk3 UjzVii — i?f 2 UjiV i2 — R kl Uj 2 Va — h T (R^ 3 Uj 2 V i3 — R kl Uj 3 V i3 ) 



The qi-q-xl couplings are given by 



n Q _ <nQ Af 



u ik ~ /x -in u kn 
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with 



f f 

JLk 
a Rk 



Bl 



h f 

a Lk 

f f 
JRk 



and 



a Lk — 
flk = 
a Rk — 
fkk = 



-h t N, 



V2 
2 



N k2 + \ g -N kl 



-htNi 



kA 



ki 



h b Lk = 


-h b N k3 


tb _ 

JLk — 




a Rk — 


-h b N k4 


fb _ 
JRk — 





where is the mixing matrix of the neutralinos. The couplings U-bj-W + read 

A w_ = ( A w Y = —( cos9 b cos9 i -sin^cos^- 
ttb] bttj V% \ — cos 9i sin 9} sin 9~ h sin 6+ 

The couplings U-bj-S^ are given by 




where 



An = v^Rlt +v 2 h 2 t R s k 2 -^g 2 J2 Vj R s k j 

A 12 = V2h b (A b Rft + fiRg - v z h T R s k t) 

A 2 i = V2ht(A t Rj£ + fiR s k t - e 3 R^) 

A22 = h b h t {v 2 R s k t +viRS) 
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The couplings t-b-S + are given by 



Ctb - t{htRkv Pl + h b R%x Pr 



The W + ~Xj ~Xk couplings read: The W + ~Xj ~Xk couplings read: 



OH 
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V2 



N M + V n N k2 



U 



J 2 



V2 



N k3 + U n N i2 . 



The S + ~Xj ~Xk couplings are given by: 



Qijk 



K2 [ V 3 iN j4 + V j2 9 -N a +N i2 )+h T (R S k tV jS N i5 + 1$; V : \ , 



+V2N tl V j3 R 



hi 



n R 1 

^ikj 



-9<RS 



U n N t3 - % ( 9 -N iX + N i2 



- h T ( R s k tu j2 N l5 + R s k tu j3 N i3 



1 

V2 



RSu j3 (^N a + N i2 
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